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Abstract 

Biskup et al. [Phys. Rev. Lett. 84 (2000) 4794] have recently suggested that the loci of partition 
function zeroes can profitably be regarded as phase boundaries in the complex temperature or field 
planes. We obtain the Fisher zeroes for Ising and Potts models on non-planar ( "thin" ) regular random 
graphs using this approach, and note that the locus of Fisher zeroes on a Bethe lattice is identical to 
the corresponding random graph. Since the number of states q appears as a parameter in the Potts 
solution the limiting locus of chromatic zeroes is also accessible. 



1 Introduction 



The idea that studying the zeroes of the partition function of a statistical mechanical model when a 
physical parameter was extended to complex values might provide valuable insights on critical behaviour 
in lattice models was first promulgated by Lee and Yang for field-driven transitions and later by 
Fisher and others for temperature-driven transitions ^ . Although there has been a considerable body 
of work studying the partition function zeroes of the Ising and Potts model on various regular lattices, 
most notably in recent years in |^ , there has been little investigation of the spin models on lattices 
exhibiting some form of geometrical disorder, such as random graphs, both for the annealed and quenched 
case. Some cases that have been investigated such as Penrose tilings have revealed intricate structure for 
the Fisher (complex temperature) zeroes away from the real axes, while still possessing Onsager scaling 
behaviour at the physical critical point Others, such as the Ising model on annealed ensembles of 
planar (j)'^ random graphs and the dual quadrangulations possess Fisher zeroes which do lie on curves in 
the complex temperature plane but which are different from the double-circle locus of the regular 
square lattice. This in itself is interesting, since it shows that a distinct locus still exists when a sum 
over a class of graphs is folded into the partition function along with a sum over spin configurations. 
For both planarfj and generic, non-planar random graphs perturbative calculations strongly suggest that 
the Lee- Yang (complex field) zeroes still lie on the unit circle in the complex activity plane 0, even 
though the technical assumptions of the original Lee- Yang theorem no longer apply in this case. 

The previous work on both fat and thin random graph models has relied entirely on perturbative 
expansions and Monte-Carlo simulations to investigate the various loci of zeroes. There is rather more 
information at hand in both these cases, however, since exact evaluations of the partition function are 
possible in the thermodynamic limit using matrix integral/orthogonal polynomial methods in the planar 
case and a straightforward saddle-point calculation in the non-planar case. In the case of the Ising model 
on regular two-dimensional lattices the particular form of the solution for the free energy, as a double 
trigonometrical integral over a logarithm, allowed a direct identification of the loci of zeroes as the curves 
or regions where the argument of the logarithm was zero. The fat and thin random graph solutions no 
longer have this form - they are even simpler, being directly the logarithms of algebraic expressions, 
which differ in the high and low temperature phases but match at the transition point, as they should. 

The general question of how to extract loci of zeroes for models with first order transitions when one 
has some analytical expression or approximation for the partition function or free energy available has 
recently been addressed by Biskup et al. They showed that under suitable technical conditions the 
partition function of a d-dimensional statistical mechanical model defined in a periodic volume V = L'^ 
which depends on a complex parameter z can be written in terms of complex functions fi{z) describing 
k different phases as 

Z = ^g,e-'3/'^ + 0(e-^/^°e-^/^) (1) 

where qi is the degeneracy of phase Z, /3 is the inverse temperature and Lq is of the order of the correlation 
length. The various // are to be interpreted as the metastable free energies of the phases, with 3?/; = / 
characterising the free energy when phase I is stable. The zeroes of the partition function are then 
determined to lie within 0{e~^^^°) of the solutions of the remarkably simple equations 

f3V{Qfi,L - 3/™,l) = TT mod 27r (2) 

where /f^'^ = fi ~ (/3^)~^ log^?;. Since we have written an expression for Z that incorporates phase 
co-existence we are implicitly looking at models with first order transitions, the canonical example being 
the field-driven transition for the Ising model. 

^We denote the planar graphs as "fat" graphs since they can be thought of as being generated by the Feynman diagrams 
in a matrix theory which has fat, or ribbon-like, propagators that carry one matrix index on each side. The fatness of the 
propagators allows the identification of a surface with a given graph and helps enforcing planarity. The converse is that 
generic random graphs appear as the Feynman diagrams in a theory with "thin" propagators that carry no edge labels and 
which have no surface interpretation. 
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The equations show that zeroes of Z asymptotically lie on the complex phase coexistence curves 
^fi,L = ^fm.L in the complex z-plane. Indeed, the idea that the loci of zeroes might be thought of as 
Stokes lines which separate different asymptotic behaviours and that the real part of the free energy in 
the various phases should match along the lines is already implicit in 

It might appear erroneous to attempt to apply these results, formulated for first order transitions 
and based on ideas arising from phase co-existence, to models where the physical temperature-driven 
transition is continuous. However, when considered in the complex temperature plane such transitions 
will be continuous only at the physical point itself and possibly some other finite set of points. This can 
be seen by looking at expressions for the magnetization for the Ising model, for example, on the square 
lattice, on fat cj)'^ graphs and on thin (p^ graphs: 



where u ^ ~ exp(— 4/3). The expressions for M apply through the complex extension of the low- 
temperature phase, and M will be zero outside this region. We can thus see that although M will vanish 
continuously at the physical critical points, u = 3 — 2\/2, c — 1/4, c— 1/3 respectively, it will generically 
be non-zero at the phase boundary approaching from within the low temperature region, whereas it will 
be zero approaching from outside. There are further discrete points where the magnetization vanishes 
continuously: at the anti- ferromagnetic point m = 3 -I- 2\/2 and the unphysical point w = — 1 on the 
square lattice; and the unphysical point c — —1/4 on the planar (f)'^ graphs. It should also be remarked 
that the apparent singularities in the expressions for Af at u = 1, c = -\/3/8 and c = 1/2,-1 respectively 
are spurious since they lie outside the domain of applicability of the expressions. Since we only have a 
discrete set of points along the phase boundary where M vanishes continuously, a typical point will be 
first order in nature and the arguments of may safely be deployed to determine the locus of zeroes. 

The results of |^ give a very natural interpretation of the loci of partition function zeroes as phase 
coexistence curves for suitably complex extended phases, both for Lee- Yang and Fisher zeroes. Such a 
picture is also in accordance with the "dynamical systems" determination of partition function zeroes 
pC| ] for models which may be defined in terms of some sort of recursive renormalization group map. In 
these cases the loci of zeroes appear as the Julia set of the map ~ in other words, the boundaries of the 
basins of attraction for the map, or phase boundaries of the complex extended phases. 

In this paper we find the locus of Fisher zeroes for the Ising model on 3-regular and 4-regular non- 
planar random graphs [jll], ^ , where such an expression can be calculated by saddle point methods and 
on the corresponding Bethe lattices where rather different recursive methods arrive at an identical answer. 
As a warm up exercise we show that the methods employed in the study of 3-regular and 4-regular graphs 
also reproduce the known results for zeroes that have been obtained by transfer matrix methods [ll3| in 
the one-dimensional |^ Ising model. The g-state Potts model on 3-regular and 4-regular random lattices 
may be treated in a similar manner to the Ising model and the locus of Fisher zeroes obtained for general 
q. Since q appears in the solution for the q state Potts model as a parameter we are also at liberty to 
look at chromatic zeroes for our random graphs, which involves looking at the saddle point solutions for 
the ground state of the anti-ferromagnetic Potts model in the complex g-plane. 

2 Ising model in one dimension 

The partition function for the Ising model with Hamiltonian 



M 



M 



M 



(l-m)i/4(i-6M + M2)V8 
(1 -u)i/2 

3(1 - 16cy/^ 
3-8c2 
(l-3c)i/2 



(3) 



(l-2c)(l + c)i/2 




(4) 



<ij> 



^i.e. on a 2-regular graph. 
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where Ui = ±1, on 2-regular random graphs with n vertices (i.e. on a one-dimensional ring) may be 
written as [|l], ID 

where Nn is the number of undecorated rings and K is defined by 

= ( I 7 ) ■ («) 

The action itself is a direct transcription of the matrix model action jl^ to scalar variables 

S - \Y.^aK-,'<^,~l{z<^l^Ul). (7) 

with z — exp(/i) and c = exp(— 2/3) . In the above (f)\_ can be thought of as representing "up" spin vertices 
and (/)^ "down" spin vertices. Since the integral is gaussian we are not obliged to pursue the saddle-point 
evaluation that is necessary in the next section for 3-regular and 4-regular graphs - we can simply carry 
out the integral directly. Scaling away the vertex counting factor of g we find that the free energy per 
site is given by the logarithm of the dominant eigenvalue of the quadratic form in the propagator, 

/ - log(A±) 

^± = 2-z • 

The partition function zeroes and hence free energy singularities appear when the square root in the 
above expression is zero and both roots can contribute 

{1- z'^)^ + 4:c^z^ = (9) 

which recovers the equation for the Lee- Yang edge derived in p^ . 

3 Ising model on thin random graphs 

The partition function for the Ising model with the Hamiltonian of equ. on 3-regular random graphs 
with 2n vertices may be written as |lll n2| 



2Tn J g^''-+^ J 27rVdetA' 



(10) 



where Nn is now the number of undecorated (no spin) (jy^ graphs. K is defined as above and the action 
now has the appropriate potential for generating (j)^ graphs 

S^]^Y.^,^ab^b~l{4>l+<l>'-^)- (11) 
a, 6 

It is necessary to include the counting factor Nn to disentangle the factorial growth of the undecorated 
graphs from any non-analyticity due to phase transitions in the decorating spins. One is also obliged to 
pick out the 2n-th order in the expansion explicitly with the contour integral over g as, unlike the planar 
graphs of two-dimensional gravity, g cannot be tuned to a critical value to cause a divergence. 

The mean-field Ising transition manifests itself in this formalism as an exchange of dominant saddle 
points. Solving the saddle point equations dS/d(f)± = 

^- = ^-+c<f>+, (12) 
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which we have rescaled to remove g and an irrelevant overall factor, we find a symmetric high-temperature 
solution 

0+ = 0- = 1 - c, (13) 

which bifurcates at c = 1/3 to the low-temperature solutions 



1 + c + Vl - 2c - 3c2 



^_ = l + ^-Vl~2c-3c^_ ^^^^ 

The bifurcation point is determined by the value of c at which the high- and low-temperature solutions 
for (j) are identical, which appears at the zero of the Hessian <let{d^S/d4)^). The magnetisation order 



parameter for the Ising model can also be transcribed directly from the matrix model 1 15 



M = -± ^ 15 

and shows a continuous transition with mean-field critical exponent {(3 = 1/2). The other critical expo- 
nents may also be calculated and take on mean-field values. 

If one carries out a saddle point evaluation of equ. (|l^), the leading term in the free energy will be 
given by an expression of the form / ~ log S", where S is the action (^l]) evaluated at either the low- or 
high-temperature saddle point solutions for </) in eqs. ( p^ , p^ ). Since we wish to match 3?/ between phases 
in order to find the loci of the partition function zeroes we are demanding in this case 

\~Sl\ = \Sh\, (16) 

where 5*^ is the low-temperature saddle-point solution and Sr is the high-temperature saddle-point 
solution, both to be taken for complex c 

sh = ^(l-c)^ 

Sl = \{\ + cf{\~2c). (17) 
6 

The resulting locus of Fisher zeroes is shown in Fig. 1, and is slightly unusual in that it is not a 
closed curve as is usually the case for Ising and Potts models on regular lattices, but rather a cusp. 
It is not unique in not being a closed curve, however, since the the locus of the Fisher zeroes for the 
one-dimensional Ising model is all of the imaginary axis ^. For one-dimensional Potts models a similar 
result holds, with the zeroes lying on a vertical line in the exp(— 2/3) plane passing through the real axis 
at 1 — g/2. The random graph Ising model has mean-field exponents because of the locally tree- like 
structure, so the scaling relation relating the angle "0 which the zeroes make as they approach the real 
axis, the critical heat exponent a and the two critical amplitudes A^, A_, 

tan 2 - a)V' = — , 18) 

sm(7ra) 

predicts^?/; = 7r/4. This is indeed what is observed at the cusp of the locus (c — 1/3) in Fig.l. 

It should be remarked that a little care is needed in applying equ. ( [l^ ) since the presence of the moduli 
signs means that one runs of risk of picking up "wrong sign" solutions. This is easily resolved when the 
complex phases are extensions of real phases as in the present case, when the behaviour on the real axis 
can be examined to make sure that the correct sign has been chosen. 

^The one-dimensional locus of Fisher zeroes can be found by applying the field/temperature duality unique to the 
one-dimensional case to the Lee- Yang zeroes on the unit circle. 

*a = 0, A+/j4_ = 0. Note that the original discussion was in terms of u = c^, but since we are dealing with conformal 
transformations any conclusions about angles will still hold, as can easily be confirmed by plotting the locus directly in 
terms of u. 
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The determination of the Fisher zeroes locus on other random regular graphs proceeds in the same 
manner as on the (j)^ graphs. On (j)'^ (four-regular) graphs, for instance the action becomes 



S=\Y. <l>aK-,'<l>, - f « + <Ai), (19) 

and the saddle-point solutions are 



2 



\/r^, (20) 



and 

</>+ 



4c 



^(V2 + 2Vl-4c2)(l- v'l-4c2), 



i(V2 + 2yr^) (21) 



at high and low temperatures, respectively. These lead to the saddle point actions 

1 1 2 

Sl = \-\c\ (22) 

which are matched in the same manner as for the <fi graphs and also lead to a cusp locus, with the cusp 
now at c — 1/2. 

Both the (j)^ and cj)"^ curves appear rather like a Folium of Descartes with the loop that would lie to 
the right of the cusp excised. The (f)^ curve turns out to have a rather more complicated formula 

27 2 _ -(1 + fx +flx2 + fa:3) + (1^3^)3/2 

16^ ~ (l + 4a;) ^^"^^ 

which is deduced from expressing \Sh \ = \Sl\ in terms of the real, x — div, and imaginary, y = 3v parts 
of w = c — 1/3. However, the expression for the curve 

where x = y = ?sv and = c — 1/2 is indeed a rescaling of the standard Folium. 

4 Ising model on Bethe lattice 

We now turn to the Ising model on a Bethe lattice, with co-ordination number m. The free energy per 
site is arrived at in an apparently very different manner to the random graphs by considering a recursion 
relation on the branches of the tree and is given by 

(3f = ~m(3 - ^mln(l - c^) + ^ \n[c^ + 1 - c(a; + a;-^)] + ^(m - 2) ln(a; + a;"^ - 2c), (25) 

where c is again equal to exp(— 2/3) and x, defined implicitly by 

^ 1-1- g-2/(-2/32.m-l ' ^20) 

emerges as the fixed point of the recursion. This model has a phase transition at /3c = | ln[m/(m — 2)] 
and the magnetisation per site is 

M= \, , (27) 
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which ranges from —1 for a; ^ cx) to +1 for x ^ 0. The derivation of these equations can be found in 

If we restrict ourselves to ft. = and m = 3 in order to compare with the Fisher zeroes on (j)^ graphs, 
the solutions of equ. ( p6|) are given by 

a; = 1, 



l-c±Vl-2c-3c2 
. - ^ (28) 

with a; = 1 being a high-temperature solution (since AI = in this case) and the others being low- 
temperature solutions when T < Tc- If we substitute these into equ. ( p5|) we find 

PfH = -^/3-^log(l-c2) + ilog[2(l-c)3], 

Ph = -^/3-^log(l-c2) + ilog[(l + c)2(l-2c)], (29) 

where we have deliberately kept the terms that do not depend on x, and hence are the same in both the 
high- and low-temperature phases, separate. It is clear from these expressions that demanding 3? /l = 'StfH 
gives us back exactly the same expression for the Fisher locus as in the random graph case, namely: 

|2(l-c)3| = |(l + c)2(l-2c)|. (30) 

The first two terms in each of eqs. (^9|) have been spirited away into the similarity sign in / ^ log S in 
the random graph expressions, where they arise from the normalising detK in equ. (|lO|). 

That the Bethe lattice with three neighbours should have the same Fisher zeroes as the (p^ random 
graph is not unexpected, since the ratio of the saddle point solutions for the random graph model gives 
precisely the recursion relation ( p6|) and both exhibit mean-field critical behaviour. The fixed point of 
the recursion variable x in the Bethe lattice solution is equal to the ratio of the (p^ and 0_ solutions in 
the random graph case. This equivalence between the random graph and Bethe lattice models is quite 
generic and also applies to the Potts models which we now move on to discuss. 



5 Potts model on thin random graphs 

A similar expression to equ. ( p^ ) can be used to define the partition function for a g-state Potts model 
where the Hamiltonian can be written 

- /3i7 = /3 ^ (5,,,,^. - 1). (31) 

<i3> 

The spins ai now take on q values, so one might expect by analogy with the Ising model {q — 2) that q 
"fields" (j) might be needed. This leads to a Potts action of the form 

S-li^f~cY^^^<l>,~lj2^f, (32) 

i—l i<j i—1 

with c — l/(exp(2/3) + q — 2). The description in terms of q (pi fields turns out to be redundant at least 
as far as the saddle-point solution goes, since explicit solution for various q reveals that the (pi only take 
on two values in the dominant low-temperature solutions and all are equal at high temperature]^ It is 
thus possible to write down an effective saddle-point action in terms of just two fields </> = (pi...q-i and 

Sl = \{q - 1) [1 - c[q - 2)] cP'-^iq- 1)4? + _ _ _ (33) 

^This is rather similar to the one-dimensional Potts model where due to the permutation symmetry (for <? > 2) only two 
of the q eigenvalues in the transfer matrix contribute in the thermodynamic limit. 
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If we solve for 0, we find the low-temperature solutions 

(j) = ^ - - 3)c- ^ V(l - 2c(7 + 2c+ c2(j2 _ gc^g + Sc^), 

= i + i (q - l)c + i V(l - 2cg + 2c + c2q2 - Gc^g + Sc^) (34) 

as well as a high-temperature solution with (j) ~ <j) — <f>Q, 

0o = l-(g-l)c, (35) 

and another branch of the solution reversing the signs in front of both square roots. These leads to a 
high-temperature saddle-point action 

SH^l{l-{q-l)c)\ (36) 
6 

and two unwieldy expressions for the saddle-point action Sl in the low-temperature regime, one coming 
from each branch of the (f) solutions, 

Sl = -ll/6c3 + l/12g-l/4cg2-l/12cV-3/2c2g2_i3/4c3g2 

+ l/4c2g3 + 9/4c2g + ll/12cV - + 17/4c3g + 3/4cg - l/2c 

± l/12(q-2)(l-2cg + 2c + c\^-6c2(7-|-5c2)3/2. (37) 

The expressions for Sl become a lot more palatable when a particular q value is inserted. Just as for the 
Ising model the locus of Fisher zeroes can now be obtained by demanding \Sh\ = |<S'l| for complex c. 

A schematic graph of the locus of zeroes for general q > 2 is shown in Fig. 2 where it can be seen that 
the physical locus of zeroes, shown in bold, crosses the real axis in the physical region of the parameter 
(0 < c < l/{q—l)) at the first-order transition point Q = (1 — (g— l)^^/'^)/(q — 2). This physical locus lies 
inside a larger metastablc locus, shown dashed, which has a cusp at the spinodal point P — \/[2q — 1). 
The loop in the Ising locus, which was excised by comparing the actions along the real axis, has split into 
the two dashed mirror inner loops lying off the axis in the Potts models. Referring back to Fig. 1 we can 
see that in the Ising case the first-order locus vanishes but the cusp remains and takes over the role of 
the physical locus. At Q we would assign a — 1 for a first order transition point so the angle of approach 
to the real axis should be 7r/2 from equ.(^8|), which is what is observed. The other crossing point O lies 
outside the physical region < c < l/(g — 1). 

Perhaps a slightly fairer comparison with the Ising model is to plot the loci in the a = exp(— 2/3) plane 
which we have also done in Fig. 3. The outer ovoid is now the genuine locus of zeroes and the metastable 
locus maps onto the tri- foliate structure that lies largely inside this. In the variable a the physical region 
is < a < 1. 



6 Potts model and thin chromatic zeroes 

One particularly interesting feature of the Potts saddle-point solution is that q appears as a parameter, 
so it is possible to obtain the Fisher locus for non-integer or even complex g-values. This means that in 
principle one can also investigate the complex g-properties of the chromatic polynomials on the random 
graphs. The chromatic polynomial P{G, q) on a given graph G of order n encodes the number of ways 
that G may be g-coloured so that no two adjacent vertices have the same colour. It is related to the 
partition function of the anti-ferromagnetic Potts model Z{G,q,P ~ —00) on G in the zero-temperature 
limit, 

Z{G,q,P^-^) = P{G,q). (38) 

We can follow Shrock and Tsai jl^ and consider a limiting function W^({G'}, q) on some class of graphs 
{G}, in our case random regular graphs, 

W{{G},q)= lim P(G, <?)!/", (39) 

n— i-oo 
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which is related to the free energy (per site) of the corresponding Potts model 



/(G,g,/3 = -oo)=logiy({G},g), (40) 

and gives us that the saddle-point action is effectively the limiting function, S{G,q,P — —oo) ~ 
W{{G},q). Since (3 — > — oo corresponds to c ^ 1/(9 ^ 2) wc can write down an action in this limit 
by making a suitable substitution for c in equ. (^) 

S = ^liq- 1)^' + - (41) 

Instead of looking for high- and low-temperature saddle-point solutions of this action, we are now in- 
terested in finding the "high g" {Sh, 0i = 4>2) and "low g" (S'l, 0i ^ 02 ) solutions, which are given 

by 



1 



P2.± 



2(2 -g) ' 
(2g-3)±75^ 

2(g-2) 



2-9' 

(42) 



respectively. We can now ape the procedure for finding the locus of Fisher zeroes to determine the locus 
of chromatic zeroes by examining the curves defined by 15^1 = \Sh\ in the complex g-plane rather than 
the complex c-plane. Although presented in a disguised form this turns out to be the Limagon of Pascal 
shown in Fig. 4, adding nicely to our collection of classical plane curves. We can see that, as expected, 
there are real zeroes at g = 1 and q = 2 indicating that the three-regular random graph cannot be one or 
two-coloured. Interestingly a zero also appears on the negative real axis, which is a somewhat unusual 
occurrence. There is no zero at q = 3 since Brook's theorem in graph theory jl^ tells us that if G is a 
connected graph with maximal co-ordination number z > 2, and if G is not complete, then the chromatic 
number % of G is z. Here we have a 3-regular graph that is not complete, so x = 3 and we can avoid 
frustration completely with 3 colours. 



7 Discussion 

We have seen that both Fisher and chromatic zero loci may be traced out on thin random graphs by 
adapting the approach of Q which was originally formulated in the context of first order phase transitions. 
We argued that even transitions that were continuous at the physical critical point would generically 
display first order properties at all but a finite set of points and that the formulae of equ.(||) could thus 
still be employed. These gave a cusp-like locus, similar to an amputated Folium of Descartes, for the Ising 
Fisher zeroes and a rather more complicated structure for the q > 2 Potts models, where the Ising-like 
cusp was still present as a metastable locus. The limiting chromatic zero locus on (p^ graphs, remarkably, 
emerged as a Limagon. The Ising model on the Bethe lattice with 3 neighbours per site returned the 
same Fisher zeroes as the (j)^ graphs by a rather different route, showing again the effective equivalence 
of the models. 

We have not discussed the case of planar random graphs here, but a similar approach can be applied 
to the solutions of and compared with the series expansions and Monte-Carlo simulations of 0. 
Similarly, we have concentrated entirely on Fisher zeroes, but it is clear following that one can also 
examine the behaviour of the Lee- Yang zeroes. 

In summary, noting that the locus of partition function zeroes can be thought of as phase boundaries 
in the complex temperature plane gives a useful way of determining the the locus of Fisher zeroes in the 
thermodynamic limit for Ising and Potts models on random graphs, as well as chromatic zeroes. 
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Figure 2: A schematic representation of the locus of Fisher zeroes for a q > 2 state Potts model, in the 
complex c = l/(exp(2/3) +q — 2) plane. The dashed locus represents the metastable branch and would not 
give rise to partition function zeroes, whereas the inner loop, drawn in bold represents the true partition 
function zeroes. The labelled points are the spinodal point P = l/{2q — 1), the "true" first-order point 
Q = {1 — {q — l)^^^^)/{q — 2) and the point O at which the loci meet (which is outside the physical 
region < c < l/(g — 1)). 
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